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ABSTRACT 

We describe a technique that for the first time measures star formation rates (SFRs) in damped Lyman 
alpha systems (DLAs) directly. We assume that massive stars form in DLAs, and that the FUV radiation 
they emit heats the gas by the grain photoelectric mechanism. We infer the heating rate by equating 
it to the cooling rate measured by the strength of CII* 1335.7 absorption. Since the heating rate is 
proportional to the product of the dust-to-gas ratio, the grain photoelectric heating efficiency, and the 
SFR per unit area, we can deduce the SFR per unit area for DLAs in which the cooling rate and the 
dust-to-gas ratio have been measured. We consider models in which the dust consists of carbonaceous 
grains and silicate grains. We present two phase models in which the cold neutral medium (CNM) and 
warm neutral medium (WNM) are in pressure equilibrium. In the CNM model the line of sight passes 
through CNM and WNM gas, while in the WNM model the line of sight passes only through WNM 
gas. Since the grain photoelectric heating efficiency is at least an order of magnitude higher in the CNM 
than in the WNM, most of the CII* absorption arises in the CNM in the CNM model. By contrast, in 
the WNM model all of the CII* absorption arises in the WNM. We use the measured CII* absorption 
lines to derive the SFR per unit area for a sample of ~ 30 DLAs in which the dust-to-gas ratio has been 
inferred from element depletion patterns. We show that the inferred SFR per unit area corresponds to 
an average over the star forming volume of the DLA rather than to local star formation along the line 
of sight. We obtain the average SFR per unit area and show that it equals 10~^'^ MQyr~-'^kpc~^ for the 
CNM solution and lO""'^'^ M0yr~^kpc~^ for the WNM solution. Interestingly, the SFR per unit area in 
the CNM solution is similar to that measured in the Milky Way ISM. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

Star formation is a key ingredient in the formation and 
evolution of galaxies. The idea that the Hubble sequence 
is actually a sequence in present-day star formation rates 
(SFRs) and past star formation histories (Roberts 1963) is 
supported by results of population synthesis models (Tins- 
ley 1980) and by the use of precise diagnostics of SFRs such 
as emission- line fluxes (Kennicutt 1983), and UV contin- 
uum luminosities. According to this interpretation, late- 
type Sc galaxies are characterized by low SFRs that are 
independent of time, while in early-type Sa galaxies an 
initially high SFR decreases steeply with time (Kennicutt 
1998). Star formation may also influence galaxy evolution 
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through feedback. That is, shocks generated by supernova 
explosions heat gas in protogalaxies, thereby delaying cool- 
ing and collapse to rotating disks. This process has been 
invoked in hierarchical cosmologies in order to prevent the 
collapse of too many baryons into low-mass dark matter 
halos, and to explain the angular momentum distribution 
of current galaxy disks (Efstathiou 2000). 

A direct measurement of star formation in high red- 
shift galaxies would provide an independent test of these 
ideas. Madau et al. (1996) took a crucial ffist step in 
this direction when they reconstructed star formation his- 
tories by measuring the comoving luminosity density of 
star- forming galaxies as a function of redshift. While the 
original "Madau diagrams" exhibited a peak in star for- 
mation at z « 1—2, more recent analyses, which are based 
on larger samples of galaxies and correct for extinction 
of emitted starlight by dust, show no evidence for such a 
peak. Rather, the SFR per unit comoving volume, p',, in- 
creases by a factor of ~10 in the redshift interval z =[0,1] 
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and either remains constant out to the highest rcdshifts 
confirmed so far z « 6 (Steidel et al. 1999) or keeps in- 
creasing to even higher redshifts (Lanzetta et ah 2002). 
However, the galaxies from which these results are derived 
arc unlikely to be the progenitors of the bulk of the cur- 
rent galaxy population. Whereas the star formation rate 
per unit area for the Milky Way Galaxy is given by V'* ~ 
4x10"'^ M© yr~^ kpc~^ (Kcnnicutt 1998), the comoving 
SFR at z ~ 3 is inferred from the Lyman Break galax- 
ies, a highly luminous population of star forming objects 
in which V'* > 1 M© yr"! kpc'^ (Pettini et al. 2001). 
The Galaxy was unlikely to be this bright in the past, 
as stellar population studies predict that at z ~ 3, i/'* ^ 
2xl0~^ Mq yr~^ kpc~^ for Sb galaxies such as the Milky 
Way. In fact only ellipticals are predicted to have ■^^ ~ 1 
Mq yr-1 kpc-2 at z 3 (see Gcnzcl et al. 2001). This is 
consistent with other independent lines of evidence such as 
strong clustering (Adelberger et al. 1998), which suggests 
the Lyman Break galaxies evolve into massive ellipticals in 
rich clusters. As a result, the published Madau diagrams 
need not constrain the star formation history of normal 
galaxies or their progenitors. 

This is the first of two papers in which we present a new 
technique for measuring SFRs in DLAs. The idea is to 
infer the SFR from the rate at which neutral gas in DLAs 
is heated. We determine the heating rate; by equating it 
to the cooling rate; i.e., we assume steady-state conditions 
(sec § 4). We estimate the cooling rate from [C II] 158 ^m 
emission, the dominant coolant for the Milky Way ISM 
(Wright et al. 1991). Specifically, we measure the cool- 
ing rate per H atom from the strength of C II* 1335.7 
absorption. As we shall show, it is plausible to assume 
the heating rate is proportional to the mean intensity of 
far UV (FUV) radiation, which is proportional to V'* for 
a plane parallel layer. The goal of this paper is to de- 
termine the mean star formation rate per unit physical 
area, < ip*{z) >i for a given redshift bin. In the second 
paper (Wolfe et al. 2003; hereafter Paper II) we combine 
< tp* (z) > with the incidence of DLAs per unit redshift to 
obtain /9'*(z) for our DLA sample. The advantage of our 
method is this: Because DLAs are not drawn from a flux 
limited sample of galaxies, we are able to derive i/i* val- 
ues far below those determined from radiation emitted by 
rapidly star forming objects such as Lyman Break galax- 
ies. Rather, we will show that our technique is sensitive 
to V* as low as lxlO~"^ Mq yr~^ kpc~^; i.e., to SFRs per 
unit area below the flux thresholds of detectors on 10 m 
class telescopes. 

A further advantage of our technique is that it probes 
the physical state of neutral gas at high redshifts. Specifi- 
cally, from our determination of the heating rate we com- 
pute the thermal equilibrium of the neutral gas. By anal- 
ogy to models for the ISM (here and throughout the pa- 
per ISM refers to the interstellar medium of the Milky 
Way Galaxy) our calculations predict a two-phase medium 
comprising a cold (T 80 K) neutral medium (CNM) and 
a warm (T ~ 8000 K) neutral medium (WNM) in pres- 
sure equilibrium with each other. We consider two models: 
one in which the line of sight to the QSO encounters gas 
in both CNM and WNM phases and another in which it 
encounters only the WNM phase. The WNM model is 



considered because of recent arguments that DLAs consist 
only of WNM gas (Norman & Spaans 1997; Liszt 2002; 
and Kanekar & Chengalur 2001). We find that the tech- 
nique is insensitive to the masses and sizes of individual 
DLAs. This has the advantage that the results are not crit- 
ically dependent on model assumptions such as the mass 
or length scale of the dark matter mass distribution (e.g. 
Prochaska & Wolfe 1997; Haehnelt et al. 1998). 

The paper is organized as follows. In § 2 we describe 
the basic strategy for inferring the SFR per unit area from 
measurements of C II* 1335.7 absorption lines. We present 
data for 33 DLAs obtained mainly with the HIRES echelle 
spectrograph (Vogt et al. 1994) on the Keck I 10 m tele- 
scope in § 3. § 4 presents two- phase models for the neu- 
tral gas in DLAs. We explain how radiative excitations 
cause the [G II] 158 /xm emission rate to exceed the 158 
/zm cooling rate. In § 5 solutions to the transfer equation 
for FUV radiation arc given. We use these solutions to 
predict heating rates as functions of tp* and dust-to-gas 
ratio. For each DLA we measure ■^^ by combining mea- 
surements of heating rate and dust-to-gas ratio with the 
solutions. We then determine <V'*> for two redshift bins 
drawn from the full sample of 33 DLAs. The significance of 
these measurements is they refer to objects representative 
of the protogalactic mass distribution; i.e., objects likely 
to evolve into normal galaxies. In § 6 we test the assump- 
tions of our dust models for self consistency. A summary 
and concluding remarks are presented in § 7. 

2. THE IDEA 

Our technique is based on the idea that massive stars 
form out of gas in DLAs. Evidence for this stems from the 
physical resemblance between DLAs and the neutral gas of 
the ISM, the presence of heavy elements in DLAs, and the 
fact that DLAs in the redshift interval z = [2,3] contain 
sufficient baryons in the form of neutral gas to account 
for all the visible stars in current spiral galaxies (Storrie- 
Lombardi & Wolfe 2000). Although stars likely form out of 
molecular rather than atomic gas, and molecules are rarely 
detected in DLAs (Lu et al. 1997; Petitjean et al. 2000), 
the presence of heavy elements argues for the formation 
of massive stars in DLAs. Such stars emit FUV radiation 
{hv = 6 — 13.6 eV) that illuminates dust grains known 
to be present in the gas (Pei & Fall 1995; Prochaska & 
Wolfe 2002; hereafter PW02). By analogy with the Milky 
Way ISM, a small fraction of the incident photon energy is 
transferred to photo-ejected electrons that heat the gas via 
Coulomb interactions with ambient electrons (e.g. Bakes 
& Tielens 1994; Wolfire et al. 1995 [hereafter W95]). In 
this case the heating rate per H atom at displacement vec- 
tor r is given by 

r^(r) = 10-24K(r)eGo(r) ergs s'^ H'^ . (1) 
In the last equation ii{r)=k{r)DLA/kMW, where k{r)DLA 
is the dust-to-gas ratio in the DLA at r, and kMw is the 
dust-to-gas ratio assumed for the present epoch ISM of the 
Galaxy (Bakes & Tielens 1994; see discussion in § 4.1). 
The incident FUV radiation field Go equals Air J, where J 
is the mean intensity integrated between 6 and 13.6 eV, 
and is in units of Habing's (1968) estimate of the local 
interstellar value (=1.6xl0^'^ ergs cm~^ s~^). The quan- 
tity e is the fraction of FUV radiation absorbed by grains 
and converted to gas heating (i.e., the heating efficiency); 
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e is also a function of G^T^/'^/ne (Bakes & Tielens 1994). 
For a plane parallel layer, J is proportional to the source 
luminosity density projected perpendicular to the plane; 
i.e., the source luminosity per unit area, which in the case 
of FUV radiation is proportional to -0* (§ 5.1). As a re- 
sult, we can deduce V* provided we know k, F^, and e. 
Because e is well determined for a wide range of physical 
conditions (Bakes & Ticlcns 1994; Wcingartncr & Draine 
2001a), this reduces to measuring the heating rate and 
K. Both of these are obtainable from HIRES spectroscopy 
since the heating rate can be inferred from C II* absorp- 
tion and K is determined from the abundance patterns and 
metallicity of the gas (sec § 4.1). Note, to derive equation 
(1) we assume e in high-z DLAs is the same as in the ISM. 
In other words we assume quantities determining e, such 
as the photoelectric cross-section, photo-electric ionization 
yield, kinetic energy partition function, and grain size dis- 
tribution (Bakes & Tielens 1994) are the same in DLAs 
and the ISM. 

We determine the heating rate by equating it to the 

cooling rate; i.c, we assume steady-state conditions. This 
is reasonable since the cooling times are w 10^—10® yrs 
which are short compared to the dynamical time scales for 
most model protogalaxies (see § 4.3). As a result we let 

r = nA , (2) 

where the total heating rate F includes other sources of 
heat in addition to F^, and n and A are the gas density 
and cooling function. In the ISM, cooling is dominated by 
[C II] 158 yum emission, i.e., Acii, with a luminosity 
L([C II]) = 5xl0^iQ (Wright et al. 1991). The 158 /im 
line results from transitions between the and 
fine-structure states in the ground 2s^2p term of C"*". Most 
of the emission from the Galaxy and other nearby spi- 
rals arises in the diffuse CNM gas rather than from star- 
forming regions in spiral arms, or photo-dissociation re- 
gions on the surfaces of molecular clouds (e.g. Madden et 
al. 1993). The last point is especially relevant for damped 
Lya systems where molecules are rarely detected (Lu et 
al. 1997; Petitjean et al. 2000). 

Pottasch et al. (1979) used the following expression to 
estimate the [C II] 158 jim emission per H atom from gas 
detected in absorption against background sources: 

h = N{C\r)hvulAui/N{m)ergss-^YL-^ , (3) 

where N{C II*) is the column density of 0+ ions in the 
state, A''(H I) is the H I column density, and Aui and 
hvui are the coefficient for spontaneous photon decay and 
energy of the transition. In fact, Ic is just 

the density— weighted average along the line of sight of the 
more fundamental quantity, /cr(r), the rate of spontaneous 
emission of energy per H atom at a given displacement 
vector r. That is 

_ j nm(s)lcr{Y(s))ds 
J nm{s)ds 

where 

new j-r) A^ihu^i 

'crirj — , {O) 

nHi(r) 

ncii* and nni are the volume densities of H I and C II*, 
and ds is the differential path length along the line of sight. 
Notice that lcr=^i^j /nm where j is the volume emissiv- 
ity appearing in the radiative transfer equation. We can 



measure Ic since iV(C II*) and 7V(H I) arc measurable: 
N{C IF) from C II* 1335.7 absorption and A' (H I) from 
damped Lya 1215.7 absorption. Note also that F = Zcr 
when cooling is dominated by [C II] 158 jim emission . 

3. THE DATA 

We have determined Ic for 33 DLAs. The results were 
obtained by measuring A'^(C II*) and A'^(H I) from accurate 
velocity profiles. In Table 1 we show Ic and other proper- 
ties to be \iscd in subsequent analyses. Column 1 gives the 
coordinate name of the background QSO, column 2 the ab- 
sorption redshift of the DLA, column 3 the H I column den- 
sity, column 4 the C II* column density, column 5 the iron 
abundance relative to solar (where [Fe/H]=logio(Fe/H) 
— logio(Fe/H)Q), column 6 gives the silicon abundance, 
[Si/H], and Ic is given in column 7. In cases where Fe 
absorption lines were not measured, we substituted proxy 
elements such as Ni, Cr, and Al. In cases where Si ab- 
sorption was not measured, we used S or Zn as proxies 
(see PW02 for a full description of these procedures) . The 
data for 30 of the entries were obtained with HIRES: 23 of 
these by our group (Prochaska et al. 2001), 6 by Lu et al. 
(1997), and one by Songaila & Cowie (2001). Data for the 
remaining 3 entries were acquired with the UVES spec- 
trograph on the VLT 8 m telescope (Ellison et al. 2001; 
Srianand et al. 2000; Dessauges-Zavadsky et al. 2001). 

Figure 1 shows six examples of C II* velocity profiles 
used to derive A^(C II*) in Table 1 along with correspond- 
ing low-ion resonance profiles. While the velocity struc- 
tures of the two profiles exhibit overall similarity, sta- 
tistically significant differences exist. These are evident 
in (1) the DLA toward Q0347-38 (panel c) where two 
strong velocity components arc detected at 1;=— 8 and -|-12 
km s~^ in Fe II 1608, while only the i'=-|-12km s~^ com- 
ponent is detected in the C II* 1335.7 profile even though 
there is sufficient signal-to-noise to detect the v=—'Skm s~^ 
component, and (2) the DLA toward Q2231— 00 (panel f) 
where absorption in C II* 1335.7 between —50 km s^^ and 
—20 km s~^ is not detected in Si II 1808. In Paper II we 
show these differences to be evidence for a multi-phase gas 
in which /cr(r) varies along the line of sight. 

In Figure 2 we plot Ic versus A^(H I) for the DLAs. There 
are 16 positive detections (red data points), 2 lower lim- 
its (95% confidence intervals plotted as blue points), and 
15 upper limits (95% confidence intervals plotted as green 
data points). One purpose of this plot is to illustrate pos- 
sible systematic effects such as correlations between Ic and 
A''(H I). No such correlation is evident in the data. How- 
ever, there is a tendency for the upper limits on Ic to occur 
at low H I column densities: 12 of the 15 upper limits on 
A^(C II*) occur at logioA^(H I) < 20.6 cm'^. This suggests 
that at least some of the null detections with large upper 
limits result from gas with C II* column densities suffi- 
ciently low that 1335.7 is undetectable rather than from 
low values of Ic- Other upper limits are caused by blending 
between C II* 1335.7 and Lya forest lines (as for the z = 
2.154 DLA toward Q2359— 02B). However, we cannot ex- 
clude the possibility that the remaining upper limits arise 
from Ic values substantially below the positive detections. 
The lower limits correspond to cases in which C II* 1335.7 
is saturated. 

The second purpose of this plot is to compare C II* 
emission rates in DLAs and in the ISM. Thus we plot 1^ 
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Table 1 
DLA PROPERTY 









log 


10W(H 1) 

9 \ 

rill " ) 


logioiV((J 11') 

{ cui - ) 


[l-'e/HJ 


[Ki/H] 




QSO° 




V. /<« 


f 






(ergs 


s-i Il-i) 


(1) 




(2) 




(:-!) 


(i) 




(<>) 






Q0019-15 


3, 


.439 


20, 


,92±0.10 


13.84±0.02 


-1.587±0.108 


-1.058±0.113 


-26, 


.61±0.10 
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2, 


.309 
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,37±0.05 
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2, 
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20, 
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< 12.78 


- 1.770±0.102 
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-27.24 


Q0201+ll=''' 


3, 


.387 


21, 


,26±0.10 


14.12±0.10 


- 1.410±0.113 


-1.250±0.150 


-26, 


.67±0.10 


Q0255+00 


3, 


.915 


21, 


,30d=0.05 


13.44±0.04 


-2.050±0.101 


-1.779±0.052 


-27 


.38±0.07 


Q0307-49'^ 


4, 


.466 


20, 


,67d=0.09 


< 13.59 


-1.960±0.220 


-1.550±0.12 


< 


-26.60 


Q0336-01 


3, 


.062 


21, 


,20d=0.10 


14.04±0.03 


-1.795±0.105 


-1.406±0.100 


-26, 


.68±0.10 


Q0347-38 


3, 


.025 


20, 


,63d=0.10 


13.47±0.03 


-1.623d=0.080 


-1.170±0.026 


-26, 


.68d=0.11 


Q0458-02 
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.039 


21, 


,65d=0.09 
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-1.767±0.102 
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-26.38 
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3, 


.017 


20, 


,48±0.10 


< 12.55 


-1.928±0.100 


-1.686±0.100 


< 


-27.45 


Q0836+11 


2, 


.465 


20, 


,58±0.10 


< 13.12 


-1.403±0.101 


-1.154±0.110 


< 


-26.98 


Q0951-049 


4, 


.203 


20, 


,40±0.10 


13.37±0.08 


< -2.591 


-2.618±0.104 


-26, 


.55±0.13 


Q0952-019 


4, 


.024 


20, 


,55±0.10 


13.55±0.02 


-1.863±0.126 
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-26, 


.52±0.10 


Q1104~18 


1, 


.661 


20, 


,80±0.10 


13.44±0.05 


-1.476±0.102 


-1.040±0.100 


-26, 


.93±0.11 


Q1108-07 


3, 


.608 


20, 


,50±0.10 


< 12.34 


-2.116±0.101 


-1.798±0.100 


< 


-27.68 


Q1202--07'' 


4, 


.383 


20, 


,60±0.14 


< 13.06 


-2.19±0.188 


-1.809±0.141 


< 


-27.06 


Q1215+33 


1, 


.999 


20.95±0.067 


< 13.17 


-1.702±0.085 


-1.481±0.072 


< 


-27.30 


Q1223+17 


2, 


.466 


21, 


,50±0.10 


< 14.01 


-1.843±0.102 


-1.593±0.100 


< 


-27.02 


Q1232+08* 


2 


.337 


20, 


,90±0.10 


14.00±0.10 


— 1.720±0.13 


— 1.220±0.15 


-26, 


.40±0.14 


Q1331+17^ 


1 


.776 


21, 


,18±0.04 


14.05±0.05 


— 2.058±0.41 


-1.452±0.041 


-26, 


.65±0.07 


Q1346-03 


3, 


.736 


20, 


,72±0.10 


12.55±0.11 


— 2.634±0.102 


-2.332±0.100 


-27, 


.69±0.15 


Q1425+609 


2 


.827 


20, 


,30±0.04 


< 13.33 


— 1.329±0.040 


> -1.034 


< 


-26.50 




A 

4 


00/1 

.ZZ4 


20, 


,80±0.10 


\ 1 /I 71 

> 14. / i 


— l.UyDitU.iio 


> — U.yzo 


> 


-25.61 


Q1759+759 


2 


.625 


20, 


,76±0.10 


12.80±0.05 


-1.184±0.008 


-0.786±0.011 


-27 


.52±0.11 


Q1946+76^ 


2 


.844 


20, 


,27±0.06 


< 12.46 


-2.528±0.061 


-2.226±0.060 


< 


-27.33 


Q2206-19 


2 


.076 


20, 


,43±0.06 


< 13.16 


-2.606±0.062 


-2.309±0.069 


< 


-26.80 


Q2231-00 


2, 


.066 


20, 


,56±0.10 


13.71±0.04 


-1.402±0.119 


-0.875±0.102 


-26, 


.38±0.11 


Q2237-06^ 


4 


.080 


20, 


,52±0.11 


< 12.53 


-2.140±0.167 


-1.870±0.112 


< 


-27.51 


Q2343+12*'9 


2 


.431 


20, 


,34±0.10 


12.77±0.05 


-1.199±0.100 


-0.540±0.101 


-27, 


.09±0.11 


Q2344+12''>'" 


2 


.538 


20, 


,36±0.10 


< 12.95 


-1.830±0.105 


-1.741±0.101 


< 


-26.93 


Q2348-14 


2 


.279 


20, 


,56±0.08 


< 13.21 


-2.238±0.0770 


-1.917±0.0776 


< 


-26.88 


Q2359-02A 


2, 


.095 


20, 


,70±0.10 


13.70±0.06 


-1.655±0.103 


-0.778±0.102 


-26 


.52±0.12 


Q2359-02B 


2, 


.154 


20, 


.30±0.10 


< 14.48 


-1.877±0.105 


-1.583±0.101 


< 


-25.34 



"Data from UCSD Data base unless otherwise noted 

'■Data from Lu ct al. (1997) 

<=Data from Ellison et al. (2001) 

■^Data from Dessauges-Zavadsky et al. (2001). 

''Data from Srianand et al. (2000). 

^N(Q, ir) from Songaila & Cowic (2001). 

^Excluded from minimal and maximal depletion models (see § 5.2.2) 
''Excluded only from minimal depletion model (see § 5.2.2) 



A^(H I) pairs derived for representative sightlines through 
the ISM which are shown as small blue stars (Pottasch et 
al. 1979; Gry et al. 1992). While Pottasch et al. (1979) 
did not report measurement errors, Gry et al. (1992) re- 
port 1-tT errors corresponding to ~ 0.3 dex in logio(Zc)- 
The large star was derived by dividing the total 158 /xm 
luminosity of the Galaxy by the H I mass of the disk (Hol- 
lenbach & Tielens 1999); the result corresponds to the 
density- weighted average of l^r (r) integrated over the disk 
of the Galaxy which we refer to as <{lc)iSM>- Compari- 
son between the DLA and ISM data demonstrates that Ic 
averaged over the DLA sample, i.e., < Zc >, is about 1/30 
times <{lc)iSM>- Because 158 /zm emission from DLAs 
has not been detected, the analogous quantity, which is 
the density- weighted average of Zcr(r) over the entire H I 
mass distribution of the DLAs, is miknown. Nevertheless 
the data covering 33 sightlines through DLAs strongly sug- 



gest the C II* cooling rate per H atom to be much lower 
in DLAs than in the ISM. The heating rates are therefore 
correspondingly lower. 

The ratio of the two heating rates is simply explained 
if the DLA gas is heated by the same mechanism that 
heats the ISM; i.e., photoejection of electrons from dust 
grains. If the mean intensities of FUV radiation are the 
same and the photoelectric efficiencies are the same, equa- 
tion 1 shows that the ratio of heating rates equals the ratio 
of dust-to-gas ratios, fcuLyi/^MW- Pettini et al. (1994) set 
kMW equal to the mean value in the Galaxy and estimate 
that koLA/kMW ~ 1/30 to 1/20 which approximates the 
DLA metallicity relative to solar and is remarkably close 
to the ratio of heating rates. Either this is a chance coin- 
cidence or it means that Go in DLAs is nearly the same as 
in the ISM and the paucity of grains accounts for the lower 
rate of heating in DLAs. Consequently we shall adopt the 
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Fig. 1. — Velocity profiles comparing C II* and selected resonance 
lines in 6 DLAs. In each case the resonance transition and C II* 
column density, 7V(C II*), are specified. 
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Fig. 2. — Inferred [C II] spontaneous emission rate per H atom, Ic, 
versus H I column density. Red data points are positive detections 
with error bars for 16 DLAs. Green data points show 2-cr upper 
Umits on Ic for 15 DLAs, and blue data points depict 2-a lower 
limits on Ic for 2 DLAs. Small blue stars depict detections toward 
H I clouds in ISM. Large blue star is spontaneous emission rate per 
H atom averaged over disk of Galaxy 



grain photoelectric heating mechanism. 

4. MULTLPHASE STRUCTURE OF THE DLA GAS 

To further evaluate the grain photoelectric heating hy- 
pothesis, we compute the thermal equilibrium temperature 
as a function of density for gas subjected to photoelectric 
grain heating. We adopt the treatment of W95 who cal- 
culate the two-phase structure of neutral gas in the ISM. 
In this calculation the gas is assumed to be mainly atomic 
and in a state of thermal and ionization equilibrium. We 
also include heating and ionization due to cosmic rays, soft 
X-rays, and the photoionization of C I by FUV radiation 
(cf. W95). Coohng is assumed to arise from fine-structure 
and metastable transitions in ions of abundant elements, 
from Lya, and from radiative recombination of electrons 
onto grains. Rather than repeat the W95 analysis here, 
we summarize the important points, emphasizing how the 
input physics for DLAs differs from that of the ISM. 

4.1. Elemental Abundances and Dust-to-Gas Ratios 

Elemental abundances affect both the heating and cool- 
ing rates in DLAs in different ways. Consider the heating 
rate. The rate of grain photoelectric heating, r^, is pro- 
portional to the dust-to-gas ratio in DLAs, k^LA, which 
in turn depends on the fraction of extant metals in grains. 
More specifically, r^; depends on the grain composition, 
i.e., on whether grains in DLAs are mainly carbonaceous, 
as in the Galaxy, or mainly silicates as in the SMC (see 
Weingartner & Draine 2001b). Recent evidence suggests 
that at z < 1, grains in DLAs are Galactic. The unam- 
biguous detection of the 2175 A graphite absorption fea- 
ture and the overall shape of the reddening curves suggest 
that a known DLA at z = 0.524 (Junkkarinen et al. 2002) 
and a DLA detected in a gravitationally lensed galaxy at z 
= 0.83 (Motta et al. 2002) contain Galactic dust. On the 
other hand, a search for the 2175 A feature in 5 DLAs with 
mean redshift w 2 (Pei et al. 1991) resulted in null detec- 
tions with upper limits on 2175 A optical depth signifi- 
cantly below predictions based on the relative reddening 
of QSOs with foreground DLAs. To account for this result, 
Pei et al. (1991) suggested that dust in DLAs resembles 
SMC dust which does not exhibit the 2175 A feature, pre- 
sumably because it is composed mainly of silicates. 

For these reasons we shall consider a "Gal" model in 
which dust in DLAs consists of carbonaceous grains and 
PAHs and an "SMC" model in which the dust in DLAs 
consists of silicate grains. The "Gal" model assumes that 
small (4—10 A) carbonaceous dust grains dominate the 
heating as they do in the Galaxy, and we shall use the heat- 
ing efficiency for Galactic regions rich in small carbona- 
ceous grains computed by Bakes & Tielens (1994). We 
infer the abundance of carbonaceous grains (per H atom) 
from the relative dust-to-gas ratio, k, where K=kDLA/kMW 
and kMW is the dust-to-gas ratio in the Galaxy. We de- 
termine K from the observed depletion of Fe in each DLA 
(see derivation in the Appendix). Our method assumes 
that the number of C atoms depleted onto dust grains per 
depleted Fe atom is the same in DLAs as in the Galaxy; 
i.e., 
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This is a reasonable assertion because Fe tracks C in metal- 
poor stars (Carretta et al. 2000). Indexing the C depletion 
to the Fe depletion is necessary because the C abundance 
typically cannot be measured in DLAs and Fe is the only 
clement for which the depletion level is known (by com- 
parison with Si). We also assume the size distribution of 
the carbonaceous grains containing the depleted C atoms 
follows that of Bakes & Tielens (1994). As a result, the 
grain photoelectric heating rate equals computed for 
pure carbonaceous grains and solar composition multiplied 
by K (equation 1). 

The "SMC" model assumes that silicate grains dominate 
the heating. The absence of a 2175A bump in the typical 
SMC dust extinction curve and deficit of 12/im emission 
indicate a lack of small carbonaceous grains (Sauvage & 
Vigroux 1991), so it is the remaining small silicate grains 
that dominate the heating. Therefore, we use the heat- 
ing efficiency calculated by Wcingartner & Draine (2001a) 
for regions that lack small carbonaceous grains. We in- 
fer the abundance of silicate grains from the dust-to-gas 
ratio determined from the observed depletion of iron anal- 
ogous to the equation for carbon shown above. This as- 
sumes that the number of depleted Si atoms per depleted 
Fe atom is the same in DLAs as in the Galaxy, making 
the SMC model a hybrid of SMC and Galactic conditions. 
Observations of [Fe/Si] in DLAs will allow us to check 
this assumption against reality. In calculating the gaseous 
carbon abundance and using [Fe/Si] to determine k, wc 
assume that C and Si are nearly undepleted. The possible 
contradiction between using the heating efficiency of car- 
bonaceous or silicate grains while assuming C and Si to be 
undepleted will be discussed in § 6. 

To estimate k we compute the fraction of Fe in grains. 
To estimate this fraction we need to determine the intrin- 
sic abundance of Fe and compare it to its gas-phase abun- 
dance. Although previous workers used Zn to estimate 
intrinsic Fe, we shall use Si because (1) the median ratio 
<[Si/Zn]>=0.03±0.05 for a sample of 12 DLAs indicates 
Si traces Zn (implying that Zn may trace a-enhanced el- 
ements rather than Fe-peak elements (PW02)), (2) there 
is only one case of Si depletion to date (Petitjcan et al. 
2002), and (3) Si abundances have been measured out to 
z = 4.5 whereas Zn is rarely measured at 2: > 3.3. In the 
Appendix we show that in the case of grains composed of 
Fe, 

K = 10[SVH]«„t ^10[Fe/Si],„, _ io[Fe/Si]gas^ _ (7) 

In the last equation the abundance ratios [X/Y]=\ogio{X/Y) 
— logio(^/l^)0, and the subscripts "int" and "gas" refer 
to intrinsic and gas phase abundance ratios. 

Since k, is sensitive to the level of depletion, we test 
two models. In the minimal depletion model wc make 
use of abundance patterns of DLAs deduced by PW02. 
Though they derived the median value <[Fc/Si]i,ij> = 
—0.3 from the abundance pattern of metal-poor DLAs 
which are not expected to have significant dust depletion, 
we shall adopt the more conservative value of —0.2. This 
is consistent with a enhancement expected from type II 
supernovae which dominate nucleosynthesis at high z. We 
then let [Fe/Si] gas equal the observed ratio and use equa- 
tion (7) to derive k for each DLA. Because the value of 
[Fe/Si]i„t is not yet well established, we consider a sec- 



ond, maximal depletion model in which we derive k by 
assuming [Fe/Si]i„t=0; i.e., we assume the observed de- 
viations of the Fe to Si ratios from the solar value are 
caused only by depletion. In both models we assume 
[Si/Hjinf = [Si/II]gas because of evidence that Si is nearly 
undepleted. In cases where only observational limits ex- 
ist on [Si/H]gas or [Fe/H]gas we substitute elements such 
as S or Zn for Si, and Ni, Cr, or Al for Fe. One could 
also consider the prescription of Pei et al. (1999) who 
let the dust-to-gas ratio equal the observed Fe abundance, 
i.e., K = lOP'^/^ls", but that turns out to be intermediate 
between the two models considered below; the minimal 
depletion model yields the smallest values of k and the 
maximal depiction model yields the largest values of k. 

The abundances of the elements C, O, Si, and Fe in 
the gas phase affect the cooling rate since transitions of 
C+, O", Si+, and Fe+ are the major coolants. In our 
models we use the following prescription to compute gas 
phase abundances of these elements in each DLA. First 
we equate the intrinsic abundance of Si to its measured 
gas-phase abundance; i.e., [Si/II]i„t=[Si/H]gas where the 
latter are listed in Table 1. Oxygen resembles Silicon as it 
is undepleted and an a-enhanced element. As a result we 
assume that [0/H]gas=[Si/H]j„t. While we assume that 
carbon is undepleted, we are aware this poses a poten- 
tial contradiction with the "Gal" dust model in which the 
depletion level of C is assumed proportional to the de- 
pletion level of Fe and we address this issue in § 6. In 
any case C is not an a-enhanced element, but rather is 
likely to trace Fe since [C/Fe]«0 in stars of all metallic- 
ities (e.g. Carretta et al. 2000). Therefore we assume 
[C/R]gas = [Fc/R],nt where [Fc/H],„t = [Fe/Si]„,i + [Si/H]„,t. 
For consistency with the other model abundances we com- 
pute [Fe/H]gas from [Fe/H]j„t rather than equate it to the 
observed DLA Fe abundance. In this case [Fe/H]gas=[Fe/H]j„t 
-l-logio[l-/vXlO~[''^^/^l*''*]. The results are summarized in 
Table 2. 

Table 2 
Gas-Phase Abundances 

Element logio(X/H)Q logio(X/H)g»s-logio(X/H)0 



He -1.00 

C -3.44 [Si/H],„t + [Fe/Si]i„t 

O -3.34 [Si/R]int 

Si -4.45 [Si/RUt 

Fe -4.45 [Si/H],„4 + [Fe/Si],„t 

+logio[l-KXlO-[^'=/-f^l-'] 



4.2. Heating and Cooling 

In this subsection we discuss the sources of heating and 
cooling in DLA gas. Specifically, we consider a gas layer 
subjected to heating by grain photoelectric emission, cos- 
mic ray ionization. X-ray ionization, and photoionization 
of C I. We also discuss cooling by important emission lines 
from abundant elements, and show how direct excitation 
by CMB photons and indirect excitation due to pump- 
ing by FUV fiuorescence photons cause the spontaneous 
emission rate to deviate from the cooling rate. Similar 
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discussions, that did not include optical pumping, were 
given by Spaans & Norman (1994) and subsequently by 
Liszt (2002). 

4.2.1. Heating 

The heating rate is given by 

r = + Tci? + Txii + Tci , (8) 

where is given by equation (1), and Tcr, T'xr, and Fci 
arc the heating rates due to cosmic rays. X-rays, and 
photoionization of C I by the FUV radiation field, Go- 
We ignore heating due to the integrated background from 
galaxies and QSOs as it is negligible compared to F^ for 
the range of SFRs considered in § 5; i.e., logio^/"* > —4.0 
Mq yr~-^ kpc~^. We compute F by adopting expressions 
and parameters used by W95 to model the ISM, but where 
appropriate extrapolate to physical conditions pertaining 
to DLAs. Thus, in the case of "Gal" dust wc compute F^ 
by adopting the Bakes & Tielens (1994) expression for the 
photoelectric efficiency, since we assume that the DLAs 
have the same relative distribution of small grains and 
PAHs as the ISM. In the case of "SMC" dust we compute 
Td by adopting the Wcingartncr & Drainc (2001a) expres- 
sion for photoelectric efficiency in the case of pure sili- 
cates, black-body FUV radiation, and selective extinction 
Rv=3.1. On the other hand there is no a priori reason why 
Go in DLAs should equal 1.7, the widely accepted value 
for the ISM (Draine 1978). As we show in § 5, Go cx -0* 
and the SFRs per unit area in DLAs need not equal the 
Milky Way rates. Moreover, the transfer of FUV radiation 
depends on the dust optical depth which should be lower 
in DLAs than in the ISM. We use the inferred optical to 
determine a self-consistent solution for ip^,, which reveals 
the SFR per unit area (see § 5.2). 

To compute Tcr we assume VcR=CcREh{E) where ex- 
pressions for QcR- the primary cosmic-ray ionization rate, 
and Eh{E), the energy deposited for each primary electron 
of energy E, are given by W95. These authors find <^cr 
= 1.8x10"^^ s~^ in the Galaxy. We scale this result to 
DLAs by assuming 

^^^ = ^-^><^°""( lO-^.^Mr;r-ikpc-^ )-"^ ' 

where we have used logioV'* = —2.4 Mq yr~^ kpc~^ for 
the disk of the Galaxy (Kennicutt 1998). 

To compute the effects of soft X-rays we use the W95 
expressions for the heating rate, Txr, and primary and 
total ionization rates, (^xr and ^xr- We again scale to 
DLAs by assuming all these quantities are proportional to 
■0*. W95 assume soft X-rays (photon energies exceeding 
0.2 keV) are emitted by thermal and non-thermal com- 
ponents. The thermal component comprises the hot (T ^ 
10** K) coronal phase of the ISM. The non-thermal compo- 
nent consists of extragalactic power-law radiation. These 
X-rays penetrate the outsides of CNM clouds, heating the 
gas to form the WNM (e.g. Heiles 2001). W95 assume 
the incident X-rays are attenuated by a WNM layer of 
gas with H I column density logioA^,« = 19 cm~^. For our 
analysis we shall assume logioiViu = 20 cm^^ instead. This 
is because low-density {n ^ 0.1 cm"*^) WNM gas cannot 
remain neutral at H I column densities logioA^^ < 20.3 
cm~^ due to the background ionizing radiation field which 



is about 100 times more intense at high redshifts than at 
z = (e.g. Prochaska & Wolfe 1996). Had we assumed 
logioA^TO =19.0 cm~^, the total column density of neu- 
tral plus ionized gas would exceed logA/'=20.0 cm~^, and 
it is the total column density which is crucial for deter- 
mining X-ray opacity. Furthermore, Vladilo et al. (2001) 
present strong arguments that gas in DLAs is mostly neu- 
tral. Thus our limit is conservative. Another reason for 
adopting the larger WNM column density is our model 
assumption that a significant fraction of the H I c;ohmin 
density in each DLA consists of WNM gas. The result of 
these model assumptions is that cosmic rays will be a more 
important source of heating and ionization than X-rays. 

Notice, that wc have assumed that Go, Fc-R; and Txr 
are all proportional to While unproven for Tcr and 
Txr, we believe this assumption is reasonable as all three 
quantities are ultimately driven by the formation rate of 
massive stars. This is because cosmic rays are thought to 
be accelerated in supernova remnants, and much of the 
soft X-ray emission is thought to arise in hot gas located 
behind supernova shocks (McKee & Ostriker 1977). 

4.2.2. Cooling 
The cooling rate (erg cm^ s~^) is given by 

A = AfS + ^MS + ^Lya+ ^GR ■ (10) 

At T<3000 K, cooling is dominated by the fine-structure 
term, Aps- The leading contributors are emission by the 
fine-structure lines [C II] 158 /um, which typically domi- 
nates at T < 300 K, and [O I] 63 /xm which is comparable 
to 158 /im emission only at T > 300 K. Following W95 we 
also include fine structure cooling from other transitions in 
0° (i.e., neutral Oxygen) and from transitions in Si+ and 
Fe"*". At T> 3000 K, the term Ams becomes important. 
This arises from excitation of metastablc transitions of C"*", 
0°, Si+, and S+. At higher temperatures, the Lya cooling 
term, A^ya starts to dominate along with Aqr, the grain 
recombination rate. We computed Aqr by adopting the 
Bakes & Tielens (1994) expression for cooling due to ra- 
diative recombinations of electrons onto PAHs and grains. 
Note, we have not included cooling by transitions in the 
neutral species C°, Fe°, Mg°, and Si" considered by W95 
as their contribution to A is negligible. 

By definition the cooling rate of the gas equals the net 
loss of thermal kinetic energy per unit time. In the case of 
line cooling this is the product of (a) the difference between 
the coUisional excitation and de-excitation rates, and (b) 
the energy of the atomic transition. The former equals 
the spontaneous emission rate, provided collisions arc the 
dominant source of excitation and de-excitation. As a re- 
sult nAcii = Icr in the ISM where Icr is the spontaneous 
emission rate per H atom of the ^^3/2 transition. 
However, this equality can break down in DLAs since ra- 
diative excitation can be important. At high z the CMB 
contributes significantly to the rate at which the 
and ^P3/2 states are populated. Moreover, for large val- 
ues of Go, these ground-term fine-structure states can be 
populated indirectly through FUV excitation of higher en- 
ergy states; i.e., through optical pumping (termed "fluo- 
rescence" by Silva & Viegas 2002). When radiative exci- 
tations are important we have 

Icr = ?lAciI + {lcr)pump + {lcr)CMB , (H) 
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where {lcr)pump and {lcr)cMB arc the spontaneous energy 
emission rates in the hmits of pure optical pumping and 
CMB excitation. In deriving the last equation we used 
the condition 1+z «hVui/kTcMB='^'i where hVui is the 
excitation energy corresponding to the ^-P3/2^ tran- 
sition in C II, and Tcmb='2'-1^ K, the current temperature 
of the CMB. We find that 



{^lcr)pump (C/H)r/u hv^i , 

CMB = 2(C/H)A„i/ii/„iexp 



huuil[k{l + z)TcMB] 



(12) 

where (C/H) is the carbon abundance, and A^i is the 
rate of spontaneous emission for the tran- 
sition. The quantity Ti^ is the net rate at which state 
I pumps state u (see Silva & Viegas 2002). We calcu- 
lated Icr using standard expressions for excitations due 
to collisions and CMB radiation. We used the Silva & 
Viegas (2002) code, POPRATIO, to compute the rate at 
which the and ^^3/2 states are populated by optical 

3/2 



pumping. We considered indirect excitation of the P-- 



and Pi/2 states through transitions to 8 higher levels. For 
consistency we used the spectral form advocated by Draine 
(1978) for the FUV radiation field, normalized such that 
Go=J 47r J^(ii//(1.6x 10~^ ergs cm~^ s^^). However, this 
procedure ignores the effects of line opacity, which can 
effectively suppress the pumping rate when the gas is op- 
tically thick to transitions such as C II 1334.5, and CII 
1036.3 (Sarazin et al. 1979; Flanncry et al. 1979; 1980). 
Because the values of Go, resonance- line optical depth, 
and collisional excitation rates are required to evaluate the 
suppression of optical pumping, wc shall re-evaluate our 
"optically thin" solutions for self-consistency in § 5.2.2. 

4.3. Phase Diagrams for DLAs 

In this subsection we compute the two-phase structure 
of neutral gas in DLAs. Wc find that if the gas is CNM, the 
SFRs per unit area deduced for local disk galaxies (Ken- 
nicutt 1998) generate [C II] spontaneous emission rates 
similar to those observed in DLAs for typical DLA metal- 
licities. If the gas is WNM, the required SFR per unit area 
must be at least a factor of 10 higher. 

We solved the equations of thermal and ionization equi- 
librium for gas at constant density with the numerical 
techniques and iterative procedures outlined in W95. We 
checked our technique by computing solutions for ISM con- 
ditions. In that case we assumed Go=1.7, [Si/H]i„t=0, 
logio«;=[Fe/H], and the same density-dependent depletion 
formulae advocated by W95. The results are in good, 
though not exact, agreement. Most importantly, the Icr 
versus n curves are in excellent agreement with the W95 
results except at logign < —0.5 cm~^ where optical pump- 
ing effects ignored by these authors causes Icr to deviate 
significantly above nAcii- 

To illustrate the behavior of two-phase media with DLA 
conditions we let [Si/H]j„t = —1.3, the mean Si abundance 
found for DLAs (PW02). We assume the "Gal" model 
and use maximum depletion to find login k = —1-5. We 
assume an ISM radiation field. Go = 1.7, and adopt the 
mean redshift of the DLA sample, z =2.8, to compute the 
CMB temperature. We compute the cosmic ray and X- 
ray heating rates from equation (9) by assuming the ISM 



star formation rate logio'!/'*=— 2.4 Mq kpc^^ yr~^. The 
resulting equilibrium curves shown in Figure 3 exhibit the 
same two-phase equilibria found by W95 for the ISM. In a 
plot of pressure, P/k, versus density, n, (sec Fig. 3a) the 
regions of thermal stability occur where d{logP) / d{logn) 
> (in the case of constant F). Thus a two-phase medium 
in which a WNM can remain in pressure equilibrium with 
a CNM can be maintained between P'"'"'/^ « 460 K cm~^ 
and p™°'^ /k K. 1750 K cm~^. An example in which P = 
(PminPmaxY^^ IS shown as tlic horizontal line connecting 
' the WNM and CNM. The intercepts with the P{n) curve 
in the WNM and CNM correspond to thermally stable 
states; a WNM with T » 7600 K and logion ss - 1 cm-^ 
in presure equilibrium with a CNM with T w 80 K and 
logiQJT. ~ +1 cm""^. Gas with densities — 0.6 < logion < 
0.0 cm~'^ is thermally unstable and evolves either to WNM 
or CNM states. Figure 3b shows the fractional ionization 
as a function of density. 

Figure 3c plots the heating rates, F (red curves), cooling 
rates, nA (green curvesand dotted blue curve in the case 
of C II), and the spontaneous emission rate l,.r (solid bhic 
curve) . It is evident that grain photoelectric heating dom- 
inates in the CNM while cosmic ray heating dominates in 
the WNM (see W95). By contrast to the ISM, cosmic rays 
dominate X-ray heating in DLAs for all densities, owing 
to the higher X-ray opacity of the H I column density as- 
sumed for DLAs. The dominant coolant in the CNM is [C 
II] 158 /zm radiation which is insensitive to density at 0.5 
< logion < 4.0 cm~^. This breaks down at logion > 4.5 
cm~^ (not shown) where C I photoionization dominates 
the heating rate and the cooling rate increases rapidly with 
density. On the other hand [C II] 158 /xm emission com- 
prises less than 10% of the cooling in the WNM. Further- 
more owing to CMB excitation and optical pumping, the 
population of the ^P3/2 state of C"*" is larger than in the 
case of collisional excitations and de-excitations alone. As 
a result, at logion < 0.0 cm~^ the spontaneous emission 
rate, Icr, will exceed the cooling through 158 /zm emission 
(sec equation 11). Notice that Icr is insensitive to density 
in the CNM where Icr^ F^ and in the WNM where Icr 
^ {h:r)pump (scc cquation 12). This reduces the uncer- 
tainty in estimating Icr in the models discussed in § 5. In 
any case the spontaneous emission rate in the CNM, logio 
Icr ~ —26.6 ergs s~^ H~^, is comparable to the mean Ic 
for the DLA sample in Figure 1. Stated differently, the 
hypothesis of grain photoelectric heating can account for 
the mean [C II] cooling rate of DLAs provided (a) heating 
occurs in CNM gas with a low dust-to-gas ratio, and (b) 
heating is driven by an FUV radiation field comparable to 
that inferred for the ISM. On the other hand the observed 
Ic could also arise in WNM gas (cf. Norman & Spaans 
1997; Kanekar and Chengalur 2001), provided Go is about 
a factor of 30 or more higher. We will address this further 
in §5. 

In our models the heating rates for DLAs were inferred 
from the cooling rates by assuming steady-state condi- 
tions. To determine whether this assumption is valid, con- 
sider the cooling time of gas in pressure equilibrium, 

{f,/2){l.l + x)kT 



tc 



nA 



(13) 



(see equation 10 in W95). We find that nA=3xlO ^ ergs 
g-i jj-i jg required for CNM gas to match the inferred 
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Fig. 3. — Two-phase diagrams for gas heated by grain photoelectric emission plus cosmic rays and soft X-rays, where the SFR per unit area, 
logioV'* =—2.4 M0 yr~^ kpc~^. Metallicity, [C/H]=— 1.5 and dust-to-gas ratio logio«;=— 1.7. Panel (a) shows pressure versus density. The 
S shape curve is indicative of a two-phase medium. Labels and horizontal line in (n, P) plane explained in text. Panel (b) shows fractional 
ionization versus density. Magenta curves in panel (c) show grain photoelectric heating (PE), cosmic-ray heating (CR), X-ray heating (XR), 
and CI photoionization heating rate versus density. Dotted blue curve is [C II] 158 cooling rate, and green curves are [O I], [Si II], Lyo! , 
and grain recombination cooling rates. Solid blue curve is [C II] 158 fira spontaneous energy emission rate. Black curve (F) is total heating 
rate. Panel (d) shows temperature versus density. 



[C II] emission rate. Since T« 50 K, we have tcooi ~ 
3x10^ yrs. If the DLA gas is in the WNIVI phase, then 
nA=3xlO~^^ ergs s""'^ H"^ is required to explain the ob- 
served Ic- In that case T « 8000 K and tcooi ~ 5x10^ 
yrs. Because these time-scales are comparable to the dy- 
namical time-scales of individual interstellar clouds, the 
assumption of thermal balance may break down on the 
spatial scales of typical interstellar clouds. However, the 
measured quantity for DLAs is Ic which is the density- 
weighted average of the C II* spontaneous emission rate, 
/cr(r), along the sightline through a typical DLA. In this 
case the relevant dynamical time scale is that of a typical 
protogalaxy which in any scenario is large compared to 
5x10^ yrs. Stated differently the fluctuations of heating 
and cooling rates integrated over the length scales of typi- 
cal DLAs average out so that the mean rates are equal. As 
a result, the assumption of thermal and ionization balance 
should be an excellent approximation for DLAs. 

On the other hand the assumption of pressure equilib- 
rium is not well established empirically. Accurate Arecibo 
21 cm measurements of H I spin temperatures in the ISM 
reveal strong evidence for CNM gas with T=25 — > 50 K, 
but no evidence for WNJVI gas with T > 7000 K. Rather, 



a significant fraction of the warm gas lies in the thermally 
unstable regime with r=500 5000 K (Heiles 2001). 
While these measurements do not rule out multi-phase 
models for the ISM, they bring to mind alternative sce- 
narios. Specifically, Vazquez-Semadeni et al. (2000; see 
also Gazol et al. 2001) compute 2D numerical simulations 
in which the dynamics of ISM clouds are dominated by 
turbulence rather than thermal instability. In this sce- 
nario the boundaries of dense CNM clouds are accretion 
shocks comprising thermally unstable gas rather than qui- 
escent contact discontinuities separating disparate phases 
at constant pressure. Moreover the unstable gas is found to 
comprise a significant fraction of the ISM mass (although, 
higher-resolution 3D numerical simulations by Kritsuk & 
Norman (2002) find the fraction to be less than 15%). But, 
as stressed by Vazquez-Semadeni et al. (2000), the cool- 
ing times in the thermally unstable gas are shorter than 
the dynamical time scales, and as a result the thermally 
unstable gas evolves quasi-statically through a sequence of 
thermal equilibrium states. For these reasons, the cooling 
curves shown in Figure 3c also apply to the "turbulence 
scenarios" . Although nAcii and Ic increase as T decreases 
from 7500 K to lOOOK, both quantities are still small com- 
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pared to the total heating rate, F. As a result, relatively 
large SFRs are required for scenarios in which C II* ab- 
sorption occurs in warm neutral gas (as discussed in Pa- 
per II) , whether or not that gas is in pressure equilibrium 
with the CNM. And it is this property which ultimately 
rules out the WNM models. We conclude that the possible 
break down of pressure equilibrium has little effect on our 
results. 

5. THE STAR FORMATION RATE PER UNIT AREA 

We now estimate the SFR per unit area for each of our 
sample of DLAs. We first solve the transfer equation for 
the mean intensity of FUV radiation corresponding to the 
K derived for each DLA and for a wide range of tp^. For 
each DLA we assign appropriate gas-phase abundances 
and then combine equation (1) with the steady state as- 
sumption of equation (2) to compute Icr- We compare the 
computed with the observed Ic to deduce '^^J^, for each 
DLA in both the CNM and WNM models. As we shall 
show, these SFRs are global in nature as they correspond 
to tjji, averaged over the entire DLA. 

5.1. Solutions to the Transfer Equation 

Assume the gas, dust and stars comprising DLAs are 

uniformly distributed throughout plane-parallel disks with 
half-width, h, and radius, R. A disk is a reasonable ap- 
proximation for DLAs because dissipative collapse of gas 
in all galaxy formation scenarios, including protogalactic 
clump formation predicted by CDM numerical simulations 
(Hachnclt ct al. 1998), occurs along a preferred axis, re- 
sulting in configurations resembling plane-parallel layers. 
Though uniformity is a highly idealized assumption, we 
shall show that the results do not differ qualitatively when 
this assumption is relaxed (see Paper II). We compute the 
mean intensity, J,^ , at the center and midplane location of 
the uniform disk and find: 
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where we have ignored scattering of photons, and ki, is the 
absorption opacity of dust at frequency v. The quantity 
is the luminosity density of the uniform disk, 6c=tan~^{R/h), 
and the extra factor of two out front comes from having 9 
run from to 7r/2. After integration we find that 



J. = 



1 (S./47r) 

2 hku 



h 



c^p{-k„R) - 



— exp(— fcjy/ix) 



(15) 



Note, we obtained equation (15) by assuming the radial 
distance to the edge of the disk equals R for 9c<0<Tr/2. 

To compute it is necessary to evaluate the quantities 
k^h and fc^i?; i.e., the dust optical depths perpendicular 
and parallel to the plane of the disk. Define the optical 
depth, T^, to be that of an average line-of-sight through 



the disk. At an average inclination angle of 45° we find 
that 

= T,/(2V2) . (16) 

To compute r^, we follow Fall & Pei (1989) who derived 
the following expression: 

= [A(A)/A(4400A)]fc DLA 

[iV(HI)/102icm-2] , (17) 

whore fcr)LA='*xfcAfw for "Gal" dust and 
koLA^KX [fcsMC X 10- [^''/■^l-''-'"^'] for "SMC" dust, andA(A) 
is the extinction at wavelength X=c/u. The [Si/HjgMC 
term is the silicon abundance of the SMC and appears 
because k is normalized with respect to Galactic dust. 
The photons responsible for photoelectric grain heating 
in DLAs have energies between £" = 6 eV and the Ly- 
man limit cutoff at 13.6 eV. At the characteristic energy 
£ = 8 eV (corresponding to A=1500 A) we find that kMW 
and ^(1500A)/A(4400A) equal 0.79 and 2.5 for the "Gal" 
model, and ksMC and A(1500A)/A(4400A) equal 0.05 and 
5.0 for the "SMC" model. We determine by assigning 
the appropriate k for each DLA, and by using the me- 
dian of our sample distribution of H I column densities, 
iV(HI)=0.48xl0^^cm~2. For each DLA we use the sample 
median rather than the measured value of 7V(HI) because 
A^(HI) along a single line-of-sight is unlikely to represent 
the H I column encountered by most of the FUV radiation; 
since << 1, the FUV radiation which heats the grains 
is transported across the entire DLA. Consequently, « 
0.01 for values of k typifying our sample. ^ 

These estimates of t^, imply the condition k,^h << 1 
in every case. For reasonable aspect ratios, R/h, we find 
that ki,R << 1 in most DLAs, but in metal-rich objects 
the condition k,yR > 1 may hold. Therefore, in the limits 
corresponding to most models, J„ takes on the following 
simple form: 



J, = - I S,/47r 



1 + \n{R/h) - k^R 
kyh « kyR « 1 
1 - 7 - \n{k^h) + 0.5k„h 



+ 0{k,Rf... ; 

0{k,hf 



. ...;fc^/i << l,ki,R » 1 



(18) 

where (=2p^h) is the luminosity density projected along 
the rotation axis of the disk, i.e., the luminosity per unit 
area, and 7 is Euler's constant. Equation (18) shows the 
mean intensity to vary linearly with the source luminos- 
ity per area, but to be weakly dependent on metallicity, 
reddening curve, and aspect ratio. To illustrate the pa- 
rameter dependence of Jy we solved equation (15) for (a) 
the "Gal" and "SMC" dust assumptions, (b) the maxi- 
mal and minimal depletion models, and (c) aspect ratios 
in the range prediced by current models of galaxy forma- 
tion. The results are shown in Figure 4 where we plot 
mean intensity versus [Si/H]. As expected, for fixed values 
of R/h, Jv increases with decreasing [Si/H] until the disk 

^We use equation (17) to determine Tv by adopting the expression 
in equation (7) for k to infer kni^A from either kMW or ksMC- 
Though the definition of these dust-to-gas ratios differs from that 
used to compute F^, the absolute values of koLA are used only for 
evaluating r^. To compute F<j, all we need is the relative quantity, 
K. Since we employ the same expression for k for both purposes, we 
are assuming that all the intrinsic properties of the grains are the 
same in DLAs and in present-day objects. 
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becomes optically thin in the radial direction. At lower 
values of [Si/H], equals a constant that increases with 
increasing R/h. At higher values of [Si/H], Ji, decreases 
with increasing [Si/H], becoming insensitive to R/h as the 
disk becomes optically thick: this is because cx Y.^/ti, 
in the optically thick limit. Figure 4 also plots results 
for a uniform sphere which, while unrealistic, provide a 
lower limit to the strength of the radiation field. For our 
standard model we shall adopt a disk with an aspect ratio 
R/h~10. From Figure 4 we see that for a measured [Si/H] 
the uncertainties m R/h and dust composition cause un- 
certainties in for the disk models which do not exceed ^ 
30%. When the spherical model is included these remain 
less than 50%. 

The grain heating rate is determined by J, i.e., Ju in- 
tegrated between photon energies £= 6 and 13.6 ev. As- 
suming the Draine (1978) UV spectrum we find that J^, 
= 10^19 ergs cm-2 s'^ sr^^ Hz^i at f = 8 eV (i.e., 1500 
A) results in 47rJ=1.6x 10"'^ ergs cm~^s~^, which equals 
Habing's (1968) estimate of the UV interstellar radiation 
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Fig. 4. 



Solutions to transfer equation given in equation (15). 
The resultant mean intensity Ji, (in units of 10~^^ ergs cm~^ s^^ 
sr~^ Hz~^), which we denote as Go, is plotted versus [Si/H] assum- 
ing Si, corresponding to a SFR, logio5/'*=— 2.4Mq yr~^ kpc^^. In 
Figure 4a we assume "Gal" dust and show results for various as- 
pect ratios, R/h, and for the maximal and minimal dust depletion 
models. Figure 4b is the same as 4a, except that "SMC" dust is 
assumed. Spherical solutions also shown. 



field. Therefore, we shall assume 
J. 



Gq 



1x10^1^ ergs cm~2 sr^^ Hz^i 




To relate Go to the rate of star formation we use the Madau 
& Pozzetti (2000) calibration. In that case 

= 8.4xlO"^^(?/5*/ Mq yr^i kpc"^) ergs cm'^ s'^ Hz" 

(20) 

By combining equations (15), (19) and (20), we can con- 
vert "0* into Gq. To see whether our technique reproduces 
the ISM radiation field, we solved for Gq assuming R= 
20 kpc, h — 0.125 kpc, logioK=0, and logioV'* = —2.4 
Mq yr~^ kpc~^. We found that Go = 1.6, which is in ex- 
cellent agreement with the Draine (1978) value of Go = 
1.7. 



Before changing topics we wish to emphasize several 
points. First, for typical metallicities the DLAs will be 
optically thin to FUV radiation in every direction. In the 
uniform disk approximation, sources at all distances con- 
tribute roughly equally to J^. As a result, the SFRs per 
unit area inferred from the C II* absorption profiles are 
representative of stars distributed throughout the entire 
DLA, not just in regions adjacent to the line of sight. This 
is in contrast to the ISM where the high dust opacity re- 
sults in SFRs with only local significance. Second, we com- 
puted Ji, for midplane points displaced from the center of 
the disk; i.e., we computed J^{r,Z) at cylindrical coordi- 
nates (r,0). When R/h=10 we found J^{r,0)/ J^{0,0) to 
slowly decrease from 1 at r=0 and to equal 0.9 at the 
median radius r = i?/V2, 0.75 at r = 0.9i?, and 0.5 
at r = 0.98i?; when R/h^lOO, J^(r, 0)/ Ji.(0,0)=0.94 at 
r=R/V2, 0.85 at r=0.9R, and 0.7 at r=0.98i?. There- 
fore, radiation fields computed from equation (15) result 
in heating rates representative of sightlines selected to 
have arbitrary impact parameters. Third, we computed 
Jj/(0, Z) at distance Z above the midplane of the uniform 
disk and then averaged the result along sightlines through 
the disk. This is a more realistic simulation of the depen- 
dence of the heating rate on mean intensity than estimat- 
ing Jjy at midplane. It is encouraging that the resulting 
mean intensities differed by less than 10% from the solu- 
tion in equation (15) for dust-to-gas ratios, —3.0 < logioK 
< 0.0. We also considered Draine's (1978) solution for 
Ju{R, Z) in which the radiation sources are confined to a 
uniform thin sheet at Z = 0. The same averaging process 
led to results in excellent agreement with equation (15) ex- 
cept when logioK > — 0.5 where the Draine expression fell 
significantly below our result. This occurs because Draine 
(1978) excluded sources within a critical radius in the disk 
location below the field point in order to avoid a singular- 
ity in his solution. Fourth, in realistic models of DLAs, V'* 
is not constant throughout a uniform disk with constant 
gas density, but rather is a function of r in a system in 
which gas density also changes with r. We shall examine 
the implications of this in Paper II. 

5.2. Determining in a DLA 

5.2.1. Technique 

The next step is to infer for each DLA from deter- 
minations of Ic and K. To do this we need an additional 
assumption about the physical state of the gas since Ic is 
not a unique function of the other two variables. Indeed, 
for fixed [Si/H], k, and ip^ the computed [C II] emission 
rate per H atom, Icr, varies with density (as shown in Fig- 
ure 3c). Therefore, to infer ■0* from observations we need 
to know the density of the gas. We address this problem 
by assuming the gas to be a two-phase medium with sta- 
ble CNM gas in pressure equilibrium with stable WNM 
gas. In that case the gas pressure P is restricted to lie 
between the local minimum and maximum of the pressure 
field; i.e., Pmin < P < Pmax (see Figure 3a). For a given 
0*, Icr iu the CNM is at least ten times larger than Icr in 
the WNM. Our first model, referred to as CNM, assumes 
that the typical DLA sightline encounters comparable 
H I column densities in the CNM and WNM. Because the 
empirical quantity Ic, is the density- weighted average of 
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^cr(r) along the line of sight (equation 4), Ic will thus be 
dominated by contributions from ^cr(r) in the CNM. As 
a result, Icr will be insensitive to density, as can be seen 
in Figure 3c which shows Icr to vary by less than 0.1 dex 
for 1 < logion < 4 cm~^; i.e., in the density range of the 
CNM. 

This insensitivity to density is a generic trait of CNM 
gas, as shown in Figure 5 where P and Icr are plotted 
against n for a grid of SFRs. The results for the Q0458— 02 
and Q1346— 03 DLAs are shown in Figures 5a. b and 5c, d 
respectively. These were chosen to compare results for a 
low-2; metal-rich DLA and a high-2: metal poor DLA. As 
4'* increases, P„u„ and Pmax increase in magnitude and 
shift to higher densities (see W95). We shall assume that 
the pressure of the two-phase medium equals the geomet- 
ric mean of Pmin and Pmax] i-C-i Peq — {,PminPmax)^^^ ■ In 

principle, Peq could assume any value between Pmin and 
Pmax- We were guided by Zcl'dovich & Pikcl'ncr (1969) 
who used stability arguments to derive unique solutions for 
Peq. Recent numerical simulations tend to support these 
conclusions and show that Peq is closer to Pmin than Pmax 
(Kritsuk & Norman 2002) in approximate agreement with 
our criterion. But given all the uncertainties, we shall not 
pursue the stability approach here. As discussed in § 3 
the intersection between Peq and the equilibrium curve 
Pin) results in two thermally stable roots: n=ncNM in 
the high-density CNM, and n—nwNM m the low-density 
WNM (see also Figure 3). The steeply rising black solid 
curves in Figures 5b, d connect the ticnm, lcr{ncNM) pairs 
determined for each tp^ of the SFR grid, where lcr{ncNM) 
represents Icr evaluated at the density, ucnm determined 
for a given The intersections between these curves 
and the observed Ic (shown as horizontal lines) determine 
ri'CNAi and V-"* for each DLA. These are dc;noted by "C" 
in Figures 5b, d as are the corresponding locations in the 
P,n plane. The unique lcr{n) curves passing through the 
intersection point "C" are shown as solid blue curves, and 
correspond to logioV^* = —1.90 and —2.75 Mq yr~^ kpc~^ 
for Q0458— 02 and Q1346— 03. In most cases the precise 
location of ucnm is unimportant as the Icr versus n curves 
are so flat in the CNM. Therefore, the SFRs we derive are 
relatively insensitive to expressions for Peq (see discussion 
in Paper II). As discussed in § 5.1 the heating rates av- 
eraged along the line of sight equal the heating rates at 
midplanc to an accuracy better than 10%. It follows from 
our steady-state assumption that the same is true for the 
cooling rates, and it is this which justifies the approxima- 
tion, lc=lcr- 

Figure 5 also illustrates why Norman & Spaans (1995) 

suggested that all high-z DLAs comprise neutral gas only 
in the WNM phase. Suppose Ic for the metal-poor DLA 
toward Q1346— 03 ([Si/II]=— 2.332) were increased to the 
lower limit placed on Ic for the more metal-rich DLA to- 
ward Q0458-02 ([Si/H]=-1.185); i.e., logio/c=-26.38 erg 
s"^ H-^ In that case the SFR impHed for the Q1346 
DLA would be logioV'* = — 1-2 Mq yr~^ kpc~^ which is 5 
times the SFR imphed by the same Ic for the Q0458 DLA. 
The pressure inferred for the Q1346— 03 DLA would in- 
crease from 3x10'^ K cm^'' to 10"'' K cm^'^ The Norman 
& Spaans (1997) models required high pressures at high 
redshifts, since they assumed metallicity to decline rapidly 
with redshift. They concluded that the gravitational fields 



generated by low-mass galaxy progenitors in CDM mod- 
els could not supply hydrostatic pressures as high as 10^ 
K cm~^ at z ~ 3. To solve this problem these authors 
concluded that all the neutral gas in high-redshift DLAs 
must be low pressure matter in which Peq < Pmin] i-c, the 
gas is a pure WNM. Though recent studies show that the 
metallicities of high-2; DLAs are not as low as assumed 
by Norman & Spaans (see PW02), there are other, in- 
dependent arguments for DLAs comprised of pure WNM 
gas. Specifically, the failure to detect 21 cm absorption in 
high-z DLAs with large H I column densities led Kanekar 
& Chengalur (2001) to invoke high spin temperatures as 
the explanation. Moreover, Liszt (2002) claimed the large 
CII/CI ratio detected in many DLAs as evidence for WNM 
gas. 

For these reasons we consider the alternative hypothe- 
sis that C II* absorption in DLAs originates in the WNM. 
More specifically we suppose that all low-ion transitions in 
DLAs, such as Lya 1215, Si II 1527, Fe II 1608, etc., arise 
in low density gas with T « 8000 K. According to the Nor- 
man & Spaans (1997) hypothesis this occurs because P < 
Pmin] i-e., CNM gas does not exist in high-z DLAs. How- 
ever, the detection of 21 cm absorption with spin tempera- 
tures, T, < 600 K in Q0458-02 (Wolfe et al. 1985) and 
< 1200 K in Q1331-M7 (Wolfe & Davis 1979; Chengalur 
& Kanekar 2001) rules out a pure WNM and is consistent 
with the presence of CNM in some cases. As a result we 
retain the two-phase liypotliesis, but assume the CNM cov- 
ering factor is so low that many sightlines miss the CNM 
phase and encounter only the WNM phase. Because we 
assume P = Peq , the density of the WNM would be given 
by the thermally stable root nwNM discussed above. In 
this case, the steeply rising dashed curves in Figure 5b, d 
connect the nwNM,lcr{nwNM) pairs, and their intersec- 
tions with the observed Ic are denoted by "W" . Proceed- 
ing by analogy with the CNM model we find that logioV'* 
=-0.90 and -1.75 Mq yr"! kpc'^ for the Q0458-02 and 
Q1346— 03 DLAs respectively. Because the Icr versus n 
curves are flat at densities below our solution for n^^M, 
the precise location of n^^M is not essential in this case 
either. Note, the shape of the Icr versus n curves are flat 
in the WNM because at low densities, the FUV radia- 
tion field dominates the population of the ^^3/2 and 
fine-structure states in CII through optical pumping. The 
increase of Icr with ip^ is due to the increase in the pump- 
ing rate caused by the increase in FUV radiation. By 
contrast, the high gas densities in the CNM cause Icr to 
equal the heating rate. So, in that case Icr increases with 
i/"* because the heating rate rises with increasing The 
effects of CMB radiative excitations should be recogniz- 
able at low where pumping is negligible, and at high 
z where the CMB intensity is high. They are evident in 
the {n, Ic) plane for the Q1346— 03 DLA as the sharp fiat 
cutoff at logio^c < — 28 erg s~^ H~^. A more detailed 
discussion of these effects is given in Paper II. 

Figures 5b,d show that the WNM solutions (points "W") 
require lower densities and higher SFRs than the CNM so- 
lutions (points "C"). This generic property of WNM ver- 
sus CNM solutions is a direct result of the lower fraction 
of total cooling due to C II* in the WNM models. There- 
fore, higher are required for lcr(ji) curves to intersect 
the observed Ic emission rate at the lower densities char- 



Wolfe, Prochaska, & Gawiser 



13 



4 - 



^2 



o 

00 






-24 



-26 - 



ta -28 





Q0458 

H — I 1 — I — I 1 1 1 1 1 1 1 — I h 



Q1346 

H 1 1 — I 1 1 1 1 1 1 1 — I 1 h 




-2 2 

log,o n (cm-3) 



-2 2 

log,j n (cm-3) 



Fig. 5. — Two-phase diagrams plotting P/k and Icr versus density. Results shown for DLAs toward Q0458— 02 (Figure 5a, b) and Q1346— 03 
(Figure 5c, d). The dot-dash curves in Figures 5a, c depict P{n) equilibrium solutions for the SFR grid logioV"* = — 4,— 3,...,0 TsAq yv~^ kpc~^. 
Dot-dashed curves in Figures 5b, d depict lcr{n) equilibrium solutions for the same SFR grid. The solid (dashed) black curves cutting across 
these solutions correspond to lcr{ncN m) ^^nd lcr{nwNM) for the SFR grid, where ncNM E'nd '^iViVM Eire thermally stable densities for 
which the pressure P={PminPmax)^^'^ ■ The intersections between these curves and the observed Zc (horizontal line) yield riQf^Mi ^WNMt 
and the corresponding ^/i* for each DLA. The intersections, denoted by "C" and "W", are also shown in the P{n) solutions. The two solid 
blue curves in Figures 5b and 5d are the unique solutions passing through the intersections. 



acterizing the WNM. In fact, the ?/'* derived for the WNM 
model are conservative lower limits. Had we used the P 
< Pmin criterion, the "W" intersections would occur at 
lower densities than indicated in Figure 5. But, there is 
even a more fundamental reason why higher values of ip^ 
are required for the WNM solutions; i.e., the absence of 
optical pumping. 

Sarazin et al. (1979) showed that the rate of optical 
pumping approaches zero, if the gas is optically thick in 
UV transitions depopulating the ground-term fine-structure 
states. Optical pumping occurs when an radiative up- 
ward transition l^k is followed by a downward radiative 
transition fc— >m to a different state, u. The net effect 
is to pump state u. At sufficiently large optical depths 
pumping ceases because the lines approach detailed bal- 
ance in which the pumping rates Pik=Pki and Puk=Pku- 
In that case the line intensities are given by the solu- 
tion to the two level atom, and the level populations of 
the I and u states are determined only by the popula- 
tion rates among the ground terms Giu and Gui- In the 
case of C II, the ground term is a fine-structure doublet 
and only the resonance transitions l^k need be optically 
thick for pumping to vanish. More specifically, Sarazin 



et al. (1979) show that pumping ceases when the quan- 
tity, ^= 2x[ln(T,fc)/7r)]i/2[i+(p,^,/G'/„)] is much less than 
1, where Tik and Pik specify the optical depth and transi- 
tion rates from the state to higher lying states, and 
Giu specifies collisional and radiative transition rates for 
the ^-Pi/2— > ^Pz/2 transitions. The only Z— >A: transitions 
relevant for pumping are C II 1036.3 and C II 1334.5, since 
the remaining resonance transitions, C II 858.1, C II 903.6, 
and C II 904.0, have wavelengths less than 912 A where 
Ji, will be negligible compared to the FUV intensities that 
we derive. In § 5.2.2 we compute ^ to check the WNM 
solutions for self-consistency. 

5.2.2. Results 

Using the above technique, we determined ■0* for a sub- 
set of the 33 DLAs in Table 1. The DLAs toward Q0951-04, 
Q1425-I-60, and Q1443-I-27 were omitted since in these ob- 
jects only limits were placed on [Si/H] or [Fe/H]. In prin- 
ciple, we could have evaluated k, [Si/H], and [Fe/H] at 
the limits we have determined. We rejected this proce- 
dure owing to the sensitivity of the inferred ■0* to n and 
to the sensitivity of the computed Icr versus n curves to 
metallicity. The DLA toward Q0952— 01 was excluded be- 



14 



STAR FORMATION IN DAMPED Lja SYSTEMS 



o 
a 



>% 

o 



^ -2 



00 
O 

— -3 



CNM 




. WNM ^ 


* 








X X 





-1 



-2 



-3 



Fig. 6. — DLA SFR per unit area versus absorption redshift deduced for dust with "Gal" composition and minimal depletion, (a) shows 
results for the CNM model, and (b) for the WNM model 



cause no measurement of [Si/H] or of a relevant proxy 
exists. We threw out the DLAs toward Q1759+75 and 
Q2343+12, because of evidence that C+/C < 1 and H°/H 
< 1; i.e., the gas causing damped Lya absorption is 
significantly ionized (Prochaska et al. 2002a; Dessauges- 
Zavadsky et al. 2003). This is in contrast to most DLAs 
where H°/H«l (Vladilo et al. 2001). As a result, the 
maximal depletion subsample comprises the remaining 27 
DLAs. For the minimal depletion sample we also ex- 
cluded the DLAs toward Q0201+11 and Q2344+12 be- 
cause the observed [Fe/Si] exceeds the assumed "nucle- 
osynthetic ceiling" value of [Fe/Si]i„t= —0.2. Consequently, 
the minimal depletion subsample comprises the remaining 
25 DLAs. The DLAs excluded from the subsamples are 
noted in Table 1. 

We determined for each DLA in these subsamples. 
Figure 6 shows the resulting ip^, plotted against redshift 
for both the CNM and WNM models in the case of min- 
imal dust-to gas ratios and "Gal" dust composition. For 
either the CNM or WNM models the results are quali- 
tatively similar for all four combinations of maximal or 
minimal dust-to-gas ratios and "Gal" or "SMC" composi- 
tion. Due to the lower fraction of cooling carried by [CII] 
158 /um emission in the WNM, the star formation rates 
are at least 10 times higher for the WNM model than for 
the CNM model. Though the positive detections exhibit 
an apparent decrease in tp^ with redshift in the interval 
2;=[3,4], no statistically significant evidence for redshift 
evolution exists. A Kendall tau test using only positive de- 
tections shows T=— 0.26, while the probability for the null 
hypothesis of no correlation, pKendaii =0.26. Since there 
is no empirical evidence that the DLAs with limits are 



physically distinct from those with measured ip^,, we shall 
assume that all systems are drawn from the same underly- 
ing population. The lowest upper limits on ip^ at z=3.608 
and z=4.080 are possible exceptions to this rule, as these 
points are outliers with ip^ systematically below the range 
of the main population. A self-consistent interpretation 
we shall explore is that the underlying population of DLAs 
consists of two-phase media with logioV"* ranging between 
—3 and —2 Mq yr^^ kpc^^. In most cases the sightlines 
pass through CNM and WNM gas. However, the sightlines 
through the two "outhers" (the DLAs toward Ql 108-07 
and Q2237— 06) pass only through WNM gas in which case 
the SFRs fall within the range of the underlying popula- 
tion (see Figure 6b). If future observations reduce the up- 
per limits on A'^(C II*) significantly, we would reject this 
hypothesis and attribute virtually all of Icr in these two 
systems to excitation by the CMB (see Paper II). 

Finally, we evaluated the pumping parameter, ^, for the 
25 DLAs in the minimal-depletion "Gal" sample. For 
C II 1036.3 and C II 1334.5 we found that the WNM solu- 
tions for ip^: resulted in f < 0.4 for the 11 cases of positive 
C II* detection and that f < 0.2 for 8 of these. We also 
found £_ > 0.025 for the DLA with a lower hmit on C II*. 
Therefore, lcr{n) is likely to be lower in the WNM, and ap- 
proach the cooling rate, nAcii (shown as the dotted blue 
curve in Figure 3c). When we recomputed "0* in the ab- 
sence of optical pumping, we found that for the WNM 
increased between 0.2 and 0.3 dex above the values shown 
in Figure 6b. While a more realistic treatment of optical 
pumping (see Flannery et al. 1979, 1980), is necessary 
to compute accurate values of ?/;* in the WNM, it is ob- 
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Table 3 
Average SFR per Unit Area 



Dust Model 






<-tp*> 




a 






Mq yr- 


kpc-^ 








CNM 




WNM 


0=2.15" 




2=3.70'^ 


0=2.15 




0=3.70 




(3.29±0.71)10^ 




(2.23±0.40)10--^ 


(3.39±0.66)10 


-■2 


(2.50±0.45)10-^ 


"Gal", min-^ 


(7.93±1.60)10- 


-3 


(4.52±1. 12)10-3 


(6.47±1.21)10- 


-2 


(4.38±0.86)10-2 


"SMC's, max 


(9.29±1.97)10- 


-3 


(4.45±1.02)10-3 


(4.94±0.99)10- 


-2 


(3.34±0.57)10-2 


"SMC", min 


(1.32±0.27)10- 


-2 


(6.39±1. 79)10-3 


(7.96±1.61)10- 


-2 


(4.84±0.99)10-2 



"Entries are SFRs per Unit Area 
''Mean rcdshift of low-z bin 
'^Mean redsliift of liigli-z bin 
'^Carbonaceous "Gal" dust 

^Maximal model where k = IQl^^/^Unt ^loP^/Silint - loP«/S'lgas^ , [Fe/Si]i„t=0 

■'^Minimal model where k = IQl^^/^hnt ^loP«/S'lint - loP^/Silgas j , [Fe/Si]i„t=-0.2 
9 Silicate "SMC" dust 



vious that the true values for V* lie somewhere between 
the values shown in Figure 6b and those computed in the 
limit of zero pumping rates. In what follows we adopt the 
"optically thin" solutions with pumping shown in Figure 
6b. We shall re-examine the effects of negligible pumping 
in Paper II where we compute the bolometric background 
radiation generated by the WNM solutions. 

5.3. The Average SFR per unit area, <'ijj^(z)> 

We now determine the average star formation rate per 

area <ip»> from our sample distribution of V'*- This is an 
important statistic since as we show in Paper II the star 
formation rate per unit comoving volume is proportional to 
<il>*>- Our goal is to determine <4>*> in as many redshift 
bins as possible. This is because we wish to determine the 
star formation history of DLAs, and because we wish to 
compare our results with comoving star formation rates 
inferred from flux-limited samples of galaxies for multiple 
redshift bins (e.g. Steidel et al. 1999). We find that 
dividing the data into more than two rcdshift bins results 
in statistical errors in <^*> that exceed the systematic 
errors. As a result we split the data set at the median 
redshift, z = 2.7, and determine <i/'*> in two redshift 
bins with redshift intervals 0i=[1.6,2.7] and Z2=[2.7,4.6]. 

While positive measurements of V* were inferred for 
many DLAs with detected Ic, upper limits on were 
set for an even larger number of DLAs with upper limits 
on Ic, and lower limits were set on ip^, for a single DLA 
with a lower limit on Ic- The numbers for the minimal 
and maximal depletion subsamples are 11 and 12 positive 
measurements, 13 and 14 upper limits, and a single lower 
limit, respectively. The presence of large numbers of lim- 
its among the data sets presents a challenge in estimating 
<'ip^>. The arithmetic mean is particularly sensitive to 
possible large values in the system for which only a lower 
limit could be measured. As discussed in § 5.2 we do not 



have empirical evidence that the systems with measured 
limits are physically distinct from those where the SFR has 
been detected. We therefore proceed from the assumption 
that the points with limits have been drawn from the same 
underlying distribution of SFRs as the detections. We use 
the detections to model this distribution empirically, as 
there is no consensus physical model for SFRs in DLAs 
at this redshift. We then treat the upper and lower limits 
as being drawn randomly from this empirical distribution 
truncated at the observed limit value. The mean value of 
the remaining probability distribution function is assigned 
to the data point, and the arithmetic mean of the full data 
set including upper and lower limits is then calculated. We 
have performed this calculation in two ways, using (1) the 
observed distribution of the detections and (2) a Gaussian 
in logio(V'*) fit to this observed distribution as our PDF 
which is then truncated by the observed limit. The second 
approach is designed to include a reasonable probability of 
high- valued outliers to which the mean is particularly sen- 
sitive, although the mean derived by this method does not 
differ strongly from the first method, and indeed neither 
mean differs strongly from the simple arithmetic mean of 
the detections alone. The uncertainty in the mean of our 
sample is then calculated using bootstrap resampling (cf. 
Efron & Tibshirani 1986). The bootstrap errors are larger 
than a nominal propagation of the errors on the individual 
detections because of significant scatter in the tp^, . 

We computed <tp*> with the two truncation approaches 
as well as assuming it to be the arithmetic mean of the 
positive detections. The results of the three techniques 
agree within the la errors indicating that <'(/'*> is a ro- 
bust statistic. As a result we henceforth assume <tp^> to 
be given by the mean of the positive detections. Table 3 
shows <?/'*> and the 1-a errors adopted for the CNM and 
WNM models in the two rcdshift bins, for the assump- 
tions of maximal and minimal depletion, and "Gal" and 
"SMC" dust composition. The errors are the quadratic 
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sums of the bootstrap errors and errors of individual de- 
tections. Obviously the final uncertainty in > is dom- 
inated by the systematic variation of the mean amongst 
the models. Averaging over the entries in Table 3 we 
find that logio<V'*>=-2.19l^;26 M3 yr^^ kpc'^ for the 
Cmi modcL and logio<V'*>=-1.32tE;;^? Mq yr"! kpc'^ 
for the WNM model. In Paper II we will reduce this 
systematic uncertainty by considering the consequences of 
cosmologically distributed DLAs with the derived SFRs. 

6. CONSTRAINTS ON DUST COMPOSITION IN DLAS 

Having determined the average SFRs per unit area we 
now examine the assumptings underlying the composition 
of dust to sec whether they arc self consistent. 

The heating rate, Td, is proportional to eGoUgrain/nH, 
the product of the photoelectric heating efficiency of the 
dust grains, the FUV radiation intensity, and the abun- 
dance of the dust grains that dominate the heating. In 
§4.1 and the Appendix we used the dust-to-gas ratio k 
to determine the abundance of grains that dominate the 
heating for both the Gal and SMC models. Our analysis 
implicitly assumes that the number of depleted C or Si 
atoms per depleted atom of Fe is the same in DLAs as in 
the Milky Way (see equation 6). We will now test this 
assumption for consistency with the lack of evidence of Si 
depletion in DLAs and with our assumption that Si and C 
are undepleted when calculating k and the gaseous carbon 
abundance. 

The depletion of Fc was determined by assuming that Si 
was undepleted. In that case the depleted ratio [Fe/SiJ^ep/ete 
=[Fe/Si]ga« — [Fe/Si]i„t where the intrinsic ratio [Fe/Si]j„t 
= in the case of maximal depletion and [Fe/Si]i„i = 
—0.2 in the case of minimal depletion. Since [Fe/SiJ^as 
typically equals —0.3, we have [Fe/Si]depiete — —0.3 for 
maximal depletion, and [Fe/Si]deplete = —0.1 for minimal 
depletion. Comparison with the abundance of Zn and S 
implies that Si cannot typically be depleted by more than 
0.1 dex (PW02). Over the range of [Fe/Si] observed, this 
implies that measuring the relative abundance of Fe versus 
a nonrefractory element with the same nucleosynthetic his- 
tory as Si such as S would generate values of [Fe/S] about 
0.1 dex lower than those observed for [Fc/Si]. This implies 
that our technique of using Si as an undepleted element 
results in an underestimate of the dust-to-gas ratio. This 
is because the ratio of k based on Si to k based on S is 
given by 



10[Si/s] 



1 _ loP'^/S'l'^^^IP^/S'l'"' 

1 - 10[Fc/S].,<,,-[Fc/Si],„t 



(21) 



where we assumed [Fe/Si]i„t — [Fe/S]j„f since Si and S are 
both a enhanced elements. Assuming [Si/S]=— 0.1 we find 
we have imderestimated the dust-to-gas ratios by factors 
between 1.5 (maximal depletion) and 2.2 (minimal deple- 
tion), offering the possibility of reducing the SFRs by a 
factor of two in both the Gal and SMC models (since our 
estimate of k, i.e., Ksi, is used to predict the dust grain 
abundance in both models). 

For the purposes of determining the gaseous carbon abun- 
dance in DLAs, we assumed that carbon was undepleted 
and set [C/H]gas=[Fe/H]i„t. In the Gal model, however, 
we are relying on carbonaceous grains to dominate the 



heating, so carbon must be somewhat depleted. Empirical 
determinations of the C abundance in DLAs are not avail- 
able since all detectable CII resonance lines are saturated 
(except in one system). In our Galaxy, it appears that 
at least half of the C atoms arc depleted at all densities 
(Meyer 1999). If this is also true in DLAs it would reduce 
the gaseous C abundance by half and alter our thermal bal- 
ance solutions such that <-(/'*> is reduced by w 0.2 dex; 
this is a mild change in our results that leads to no quali- 
tative differences in our conclusions. It seems more likely 
that the overall carbon depletion level is lower in DLAs 
due to the reduced metallicity. In particular, the SMC 
model has a reduction in the abundance of carbonaceous 
grains and therefore it is likely that the vast majority of 
carbon is gaseous. 

Without much knowledge of the composition of dust 
outside our own Galaxy, it is difficult to estimate the sys- 
tematic uncertainty introduced by our assumption that 
the number of depleted Si and C atoms per depleted Fe 
atom is the same in DLAs as in the Milky Way. If there 
really is a base level of C depletion independent of density 
and metallicity, this implies that using k underestimates 
the true number of depleted carbon atoms. However, this 
may still give a reasonable estimate of the small carbona- 
ceous grains which dominate the heating but do not appear 
to be part of the base depletion in our Galaxy (Sauvagc & 
Vigroux 1991). The uncertainty in the abundance of small 
carbonaceous grains is bracketed by the range of models 
for the size distribution of dust grains in Weingartner & 
Draine (2001b). The fraction of depleted C atoms in small 
carbonaceous grains could be a factor of two lower than 
implied by the extrapolated MRN size distribution used by 
our adopted Bakes & Tielens (1994) model, which would 
reduce the photoelectric heating efficiency by a factor of 
two and thereby increase our inferred star formation rates 
by a factor of two. Reducing the overall number of de- 
pleted C atoms could make the SFRs even higher, but re- 
ducing the small carbonaceous grain population by more 
than a factor of three makes silicate grains dominate the 
heating, in which case the Gal model becomes the SMC 
model. Alternatively, the fraction of depleted C atoms 
contained in small grains could be increased by up to a 
factor of four. If this is the case, or if indexing k to Fe has 
underestimated the carbon depletion, small carbonaceous 
grains could be more abundant than our assumptions im- 
ply, leading to higher heating efficiency, reduced star for- 
mation rates, and a stronger 2175A bump. Observational 
limits on the strength of the 2175A bump in DLAs make 
it difficult for small carbonaceous grains to be more than a 
factor of a few more abundant than we have supposed (Pei 
et al. 1991). Tighter observational limits on (and possibly 
detection of) the 2175A bump in DLAs arc of the utmost 
importance in reducing the systematic uncertainties in the 
nature of dust at these redshifts. 

In the SMC model, we suppose the complete absence of 
small carbonaceous grains as inferred from the lack of the 
2175A bump. This allows one to lower the depletion level 
of Si considerably and still have silicate grains dominate 
the heating. Since Fe is almost completely depleted in dif- 
fuse regions of the Milky Way, we are likely to overestimate 
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pendix would be an underestimate. This appears to be the 
case for the SMC, where absorption by clouds along the 
hnes-of-sight to Sk 108 and Sk 155 show [Si/Zn]=0, imply- 
ing a lack of silicon depletion, but [Si/Fe]=0.5, indicating 
that iron is significantly depleted (Welty et al. 1997; 2001). 
In the diffuse regions of the Galaxy modeled by Weingart- 
ner & Draine (2001b), 75% of the Si atoms and 95% of the 
Fe atoms are depleted, but in the SMC it appears that no 
more than 20% of the Si atoms are depleted even though 
70% of the Fe atoms arc;. This implies that indexing k to 
the Fe depletion overestimates the number of depleted Si 
atoms by at least a factor of three if DLAs are like the 
SMC. This uncertainty is somewhat balanced by the un- 
certainty in the fraction of Si atoms contained in small 
(< 15 A ) grains. In the models of Weingartner & Draine 
(2001b), this varies from the fraction we have assumed 
up to a factor of three higher. Therefore, the abundance 
of small silicate grains dominating the heating is unlikely 
to be more than a factor of three higher than we have 
assumed; if it were, the DLA star formation rates for the 
SMC model would be less than or equal to our current Gal 
model results. It is possible, on the other hand, to decrease 
the amount of Si depletion by an arbitrary amount, which 
would increase the inferred DLA star formation rates. For 
the WNM solution, this would worsen the conflict with ob- 
servational limits on the integrated background light dis- 
cussed in Paper II. For both WNM and CNM solutions, 
it exacerbates the general problem of overproduction of 
metals in DLAs discussed in Paper II. For the CNM so- 
lution, doubling the SMC model star formation rates for 
DLAs at z > 1.6 begins to exceed the observational limit 
on the background light as discussed in Paper II. Hence 
the integrated background light constrains the abundance 
of small silicate grains in DLAs to be at least half of that 
of the Milky Way. 

7. SUMMARY 

The conventional view of DLAs is they are high-^; neu- 
tral gas layers with low metallicities, low dust content, and 
quiescent velocity fields. In this paper we have developed 
a new technique providing a more complete picture. Using 
the C II* absorption method we find that rather than being 
passive objects transmitting light from background QSOs, 
DLAs are the sites of active star formation, and that neu- 
tral gas in DLAs is likely to be a two-phase medium. At 
this stage of our analysis it is unclear whether C II* ab- 
sorption arises in the CNM or WNM phase. Our results 
are as follows. 

(1) Our technique assumes that massive stars forming 
out of neutral gas in DLAs emit FUV radiation that heats 
the gas by ejecting photoelectrons from dust grains known 
to be present in the gas. We can infer the heating rate since 
in steady state conditions it equals the cooling rate that 
is directly measurable. This is because cooling is domi- 
nated by [C II] 158 nm emission if the gas is a cold neutral 
medium (CNM), and [C II] 158 emission per H atom can 
be obtained by measuring C II* 1335.7 absorption arising 
from the ^P3/2 excited fine-structure state in the ground 
term in C+. The heating rate equals the product of the 
dust-to-gas ratio, the mean intensity of FUV emission. Go, 
and the grain photoelectric heating efficiency, e. We can 
measure Go since the cooling rate is inferred directly from 



the column density N{C II*), the dust-to-gas ratio can be 
computed from element abundance patterns, and e is well 
determined provided the gas is CNM. By measuring Go, 
we measure the SFR per unit area. since oc Gq in 
a plane parallel layer. 

(2) We have measurements of C II* 1335.7 absorption 
in 33 DLAs; 16 of these are positive detections, 15 are 
upper limits, and 2 are lower limits. We use these data to 
infer the spontaneous energy emission rate per H atom, Ic, 
from the ratio iV(C ir)/iV(HI). We find that h in DLAs is 
typically about 1/30 times Ic measured for the ISM of the 
Galaxy. Because Ic equals the cooling rate in the CNM, in 
our model Ic equals the heating rate, which is proportional 
to KeGo, where k is the dust-to-gas ratio in DLAs relative 
to the ISM. Since k is also about 1/30, the implication is 
that Go in DLAs is similar to that in the ISM. In other 
words, V* in DLAs is similar to the local V* in the ISM, 
provided C II* absorption arises in a CNM. 

(3) We compute thermal equilibria of gas subjected to 
grain photoelectric heating and standard cooling processes. 
Since the dust content of DLAs is not well determined, we 
consider a "Gal" model in which the grains are mainly 
carbonaceous, and heating is dominated by small (< 15 
A) grains and PAHs, and an SMC model in which heat- 
ing is dominated by small silicate grains. We compute k 
from the observed depletion of Fe in each DLA and assume 
that the number of C or Si atoms depleted onto grains per 
depleted Fe atom is the same in DLAs as in the Galaxy. 
We also consider models with minimal and maximal de- 
pletion to account for the uncertainties in the Fe depletion 
levels. We include heating and ionization due to cosmic 
rays and soft X-rays. When computing cooling rates we 
account for excitation of the C+ fine-structure levels by 
CMB radiation. We also account for excitation due to 
optical pumping by FUV radiation, but find that optical 
pumping may not be significant owing to the high opacity 
in the C II resonance lines. Equating heating and cooling 
rates, we find the resulting equilibrium curves of pressure 
versus density to exhibit a maximum pressure, Pmax, and 
minimum pressure, Pmin- A two-phase medium in which 
a dilute warm neutral medium (WNM) is in pressure equi- 
librium with a dense CNM is possible if the equilibrium 
pressure, Peq, satisfies the constraint Pmin < Peq < Pmax- 
Therefore, C II* absorption can occur in the WNM or the 
CNM. Because Ic equals the cooling rate in the CNM, but 
is a small fraction of the cooling rate in the WNM, the 
SFRs implied from a measured Ic are much higher in the 
WNM than the CNM. 

(4) We calculate tp^, for selected subsets of our DLA sam- 
ple corresponding to "Gal" or "SMC" dust, and to minimal 
or maximal depletion. We first solve the transfer equation 
for sources of FUV (ss 1500 A) radiation (OB stars) uni- 
formly distributed throughout a plane parallel disk. Using 
standard reddening curves for "Gal" and "SMC" dust we 
find the disks to be optically thin parallel to the plane 
for most of our sample DLAs. As a result, the inferred 
from C II* absorption are representative of the entire DLA 
rather than just regions along the QSO sightline. To infer 
'^p^, from measurements of Ic and k we assume the equi- 
librium pressure. Peq = iPminPmaxY^'^ ■ As CXpCCtcd, WC 

find two solutions; one in which C II* absorption occurs 
in the WNM, and the other in which C II* absorption 
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occurs in the CNM. For the CNM solution we find -3.0 

< logioV'* < - 2.0 Mq yr"i kpc^^, and -2.0 < logio^, 

< - 1.0 Mq yr-^ kpc~2 for the WNM solution. Neither 
case shows evidence for redshift evolution in the interval 
z = [1.6,4.5]. In Paper II wc discriminate between these 
models by deriving cosmological constraints such as the 
bolometric background radiation. 

(5) Our assumptions that C and Si are undepletcd in 
determining the gaseous C abundance and the dust-to-gas 
ratio arc reasonable and do not create serious contradic- 
tions with the heating being dominated by carbonaceous 
or silicate grains, because the C and Si depletion levels in 
DLAs appear to be quite low. Varying the size distribu- 
tion of carbonaceous grains and the number of depleted C 
atoms per depleted Fe atom could make the SFRs reported 
for the Gal model as high as those of the SMC model, but 
one cannot lower the SFRs too far since the 2175A bump 
has not been observed in high redshift DLAs. Varying the 
number of depleted Si atoms per depleted Fe atom could 
make the SFRs reported for the SMC model as high as 
allowed by the integrated background limits, and varying 
the size distribution of silicate grains could make the SFRs 
a factor of three lower. 
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APPENDIX 

DUST TO GAS RATIO 

Consider a box containing Np°* iron atoms. If 7Vp^°* 
atoms are in the gas phase and A^]£"** atoms are locked 
up in grains, then 

Because we cannot measure iVp°* due to depletion of Fe 
atoms onto grains, we shall use Si as a proxy as Si is es- 
sentially undepleted in DLAs (see PW02). Therefore, we 
assume Np°*=Nsi{Fe/Si)int where Nsi, the total number 
of Si atoms also equals the number of gas- phase Si atoms, 
and (Fe/Si)i„f is the intrinsic (i.e., undepleted) ratio of Si 
to Fe. Let the dust-to-gas ratio k =A^^"**/7Vh- As a result 



Factoring Fe/H=(Fe/Si)(Si/H) we find 

(Si/H) (Fe/Si)i,t (Si/H)(Fe/Si)G- 



k 



^Si(Fe/Si), 



AjGas 



where A^h is the number of H atoms. Since k= k/kMw-, 
where kMW is the dust-to-gas ratio of the current Milky 
Way Galaxy, and assuming that fcMH' = (Fe/H)0, since Fe 
is almost entirely depleted, we find 

(Si/H)(Fe/Si).„t (Fe/H)g°« ^ ^^^^ 



(Si/H)0 (Fe/Si)0 
As a result 



(Si/H)0(Fe/Si)0 



10[SVH]i„ 



gas 



(A4) 
(A5) 



where we assume [Si/H]=[Si/H]j„t. 



(Fe/H)o 



(Fe/H)0 
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